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Abstract 

Q\ \ The logarithmically enhanced a 7 m(l/a) contribution to the hyperfine splitting 

of the positronium ground-state energy levels is calculated in the framework of di- 
mensionally regularized nonrelativistic quantum electrodynamics. The correction is 
negative and amounts to about 1/3 of the leading logarithmic a 7 ln 2 (l/a) one. The 
discrepancy between the experimental measurements and the theoretical prediction 
-»^. ■ is reduced. 
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Positronium, which is an electromagnetic bound state of the electron e~ and the 
positron e + , is the lightest known atom. Thanks to the smallness of the electron mass 
m e , strong- and weak-interaction effects are negligible, and its properties can be calculated 
perturbatively in quantum electrodynamics (QED), as an expansion in Sommerfeld's fine- 
structure constant a, with very high precision, only limited by the complexity of the 
calculations. Positronium is thus a unique laboratory for testing the QED theory of weakly 
bound systems. However, the theoretical analysis is complicated due to annihilation and 
recoil effects. 

The positronium hyperfine splitting (HFS) Av = E (l 3 Si) — E (l 1 ^), where E (l 1 ^) 
and E (l 3 Si) are the energy levels of para- and orthopositronium, respectively, is the most 
precisely measured quantity in positronium spectroscopy as far as the absolute precision 
is concerned. The most recent measurements of the HFS []1],0] yielded 

Az/ cxp = 203.387 5(16) GHz, (1) 
Av cxp = 203.389 10(74) GHz, (2) 

respectively. 

The present theoretical knowledge may be summarized as: 
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where Az/Q h = 7m e a 4 /12 is the leading-order result 0. The first-order correction was cal- 
culated in Ref. ||]. The logarithmically enhanced a 6 ln(l/a) term was found in Ref. f||6|]. 
The nonlogarithmic 0(a 6 ) term includes the contribution due to the radiative correction 
to the Breit potential 0, the three-, two- and one-photon annihilation contributions ||, 
the non-annihilation radiative recoil contribution ||, and the pure recoil correction com- 
puted numerically in Ref. (1(J and analytically in Ref. |H |. In 0(a 7 ), only the leading 
double-logarithmic a 7 ln 2 (l/o!) term is available [ |12]] . 
Including all the terms known so far, we have 

Av th = 203.392 01GHz, (4) 

which exceeds Eqs. (^ and (Q) by approximately 2.8 and 3.9 experimental standard 
deviations, respectively. In contrast to the well known orthopositronium lifetime puzzle,^ 
the experimental situation for the HFS is unambiguous. In fact, the experimental error 

1 For the most recent developments of this problem, see, for example, Ref. Jl3| and the references cited 
therein. 
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is compatible with a naive estimate of the theoretical uncertainty due to as-yet unknown 
higher-order corrections. Should this discrepancy persist after the dominant terms of the 
latter have been calculated, this would provide a signal for new physics. This makes the 
HFS to be one of the most interesting topics in positronium spectroscopy, both from the 
experimental and theoretical points of view. 

Thus, it is an urgent matter to improve the prediction of the HFS as much as possible, 
and one is faced with the task of analyzing the third-order correction, which is extremely 
difficult. However, there is a special subclass of the 0(a 7 ) contributions which can be 
analyzed separately, namely those which are enhanced by powers of ln(l/a) ~ 5. They 
may reasonably be expected to provide an essential part of the full 0(a 7 ) contributions. 
This may be substantiated by considering Eq. ([3D in 0(a 6 ), where the logarithmic term 
is approximately 2.6 times larger than the nonlogarithmic one. While the leading double- 
logarithmic 0(a 7 ) contribution to Eq. (||) is known [12], the sub leading single- logarithmic 
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one is yet to be found. In fact, from the positronium lifetime calculation [|Hj|.p~I 
know that the subleading terms can be as important as the leading ones. The purpose of 
this Letter is complete our knowledge of the logarithmically enhanced terms of 0(a 7 ) by 
providing the coefficient C in analytic form. 

The origin of the logarithmic corrections is the presence of several scales in the bound- 
state problem. The dynamics of the nonrelativistic e + e~ pair near threshold involves four 
different scales: (i) the hard scale (energy and momentum scale like m e ); (ii) the soft scale 
(energy and momentum scale like /3m e ); (iii) the potential scale (energy scales like f3 2 m e , 
while momentum scales like /3m e ); and (iv) the ultrasoft scale (energy and momentum 
scale like (3 2 m e ). Here j3 denotes the electron velocity in the center-of-mass frame. The 
logarithmic integration over a loop momentum between different scales yields a power of 
ln(l//3). Since positronium is approximately a Coulomb system, we have /3 oc a. This ex- 
plains the appearance of powers of ln(l/a) in Eq. ([!]). The leading logarithmic corrections 
may be obtained straightforwardly by identifying the regions of logarithmic integration 
PJT2]] . The calculation of the subleading logarithms is much more involved because certain 
loop integrations must be performed exactly beyond the logarithmic accuracy. 

In the following, we briefly outline the main features of our method developed pre- 



viously in Ref. ||13|| , where it was applied to the analysis of the subleading logarithmic 
third-order corrections to the positronium ground-state decay rates. This approach is 
similar to the one adopted in Ref. |L5 |. We work in nonrelativistic QED (NRQED) fl6| . 



which is the effective field theory that emerges by expanding the QED Lagrangian in j3 
and integrating out the hard modes. If we also integrate out the soft modes and the 
potential photons, we arrive at the effective theory of potential NRQED [[17]], which con- 
tains potential electrons and ultrasoft photons as active particles. Thus, the dynamics of 
the nonrelativistic e + e~ pair is governed by the effective Schrodinger equation and by its 
multipole interaction with the ultrasoft photons. The corrections from harder scales are 
contained in the higher-dimensional operators of the nonrelativistic Hamiltonian, corre- 
sponding to an expansion in j3, and in the Wilson coefficients, which are expanded in a. In 
the process of scale separation, spurious infrared and ultraviolet divergences arise, which 
endow the operators in the nonrelativistic Hamiltonian with anomalous dimensions. In 
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fact, these divergences completely determine the logarithmic corrections [131,15 . We use 
dimensional regularization, with d = 4 — 2e space-time dimensions, to handle these diver- 
gences, which are of the form l/e n (n = 1, 2, . . .) as e — > [11,1£,2C]. Compared to the 
traditional NRQED approach, endowed with an explicit momentum cutoff and a fictitious 
photon mass to regulate the ultraviolet and infrared behavior P, |I4| , |2~T| , this scheme has 
the advantage that contributions from different scales are matched automatically. 

In the effective theory, the HFS is generated by spin-flip operators of the effective 
nonrelativistic Hamiltonian averaged over the bound-state wave function. The hard-scale 
corrections, which require fully relativistic QED calculations and are most difficult to 
find, do not depend on (3 and do not lead to logarithmic contributions by themselves. 
However, they can interfere with the logarithmic corrections from the softer scales. The 
only results from relativistic perturbation theory that enter our analysis are the one-loop 
renormalizations of the relevant operators in the effective nonrelativistic Hamiltonian. 
The missing ingredients can all be obtained in the nonrelativistic approximation. The 
leading-order 0(a 4 ) HFS is generated by the 0{f3 2 ) spin-flip part of the tree-level Breit 
potential, which consists of the Fermi operator Vp and the annihilation operator V^ nn . The 
a 7 ln(l/o;) contribution, which we are interested in, arises from several sources. A part of 



it can be extracted from the positronium lifetime calculation [13]. This part corresponds 
to: (1) the second-order corrections in nonrelativistic Rayleigh-Schrodinger perturbation 
theory, which arise from (la) the insertions of the tree- level 0{f3 2 ) spin-independent Breit 
potential and the one-loop hard corrections to the leading-order spin-flip operators and 
(lb) the insertions of the one-loop 0{a(3 2 ) operators related to the hard [22 1, soft [19], or 
ultrasoft ]I3"| , P3"| scales (cf. Eqs. (13)-(15) of [13] , respectively) and the leading-order spin- 
flip operators; (2) irreducible corrections to the leading-order spin-flip operators, which 
include (2a) 0(a(3 2 ) and (2b) 0(a 2 (3) terms and can be obtained from the correspond- 
ing equations of Ref. ]13[ by replacing the leading-order decay operator V^p, p',S) by 
Vp + V ann . The hard renormalization coefficients of the leading-order operators Vp and 
Vann are 1 and —(44/9 + 2 In 2), respectively [pl| , p4][ . They replace the renormalization 
coefficients A P)Q of V^(p, p', S) in Ref. [|H|. Another nontrivial difference with respect to 
the positronium lifetime calculation pi is that the Fermi operator can cause a AL = 2 
transition of the spin-triplet state, so that D-wave intermediate states contribute to sec- 
ond order in nonrelativistic perturbation theory, by double insertion of Vp. This gives 
an additional contribution of 5/42 to the coefficient C. Finally, there is a modification 
of the soft part of contribution (lb) due to the spin-dependent part of the transverse- 
photon-exchange contribution, which vanishes for the annihilation channel. This gives 



an additional contribution of 68/63 to the coefficient C. Using the results of Ref. [13 
and taking these modifications into account, we obtain the following contribution to the 
coefficient C from the sources enumerated above 
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where the various contributions are given separately. Note that contribution (lb) also 
includes the entire double-logarithmic term not presented in Eq. (|5|). The structure of 
the overlapping divergences resulting in the double-logarithmic contribution to the HFS 
is similar to the positronium lifetime analysis JT3[ . 

Another part of the a 7 ln(l/a) contribution is produced in the first order of the non- 
relativistic Rayleigh-Schrodinger perturbation theory by the one-loop hard corrections to 
the 0(/3 A ) and 0(a/3 3 ) spin-flip operators, which give rise to the a 6 ln(l/o;) contribution 
to the HFS. The relevant operators are generated by the relativistic correction to the 
Coulomb-photon exchange, the relativistic correction to the transverse-photon exchange, 
the kinematical retardation, and the one-loop correction involving seagull vertex diagrams 
with one Coulomb and one transverse photon. For our calculation, we only need the three 
hard renormalization coefficients, usually denoted as cp, cs, and c PP >, which are related to 
the anomalous magnetic moment of the electron and are finite and scheme independent. 
By gauge invariance, the coefficients cp and cs are the same as for the 0{f3 2 ) opera- 
tors, while the coefficient c pp i parameterizes the new 0(af3 4 ) operator contribution to the 



transverse-photon exchange [21]. The resulting contributions to the coefficient C are 
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The remaining part of the a 7 ln(l/a) contribution also corresponds to the first order of the 
nonrelativistic Rayleigh-Schrodinger perturbation theory and is related to the relativistic 
corrections to the operators contributing to the HFS in lower orders. The relevant oper- 
ators are generated by the diagrams with one or two transverse-photon exchanges, where 
the momentum of the photon with spin- independent interaction can be soft or ultrasoft, 
and by the soft diagrams with one or two seagull vertices involving either one Coulomb 
and one transverse photon or two transverse photons. The resulting contribution to the 
coefficient C reads 

41 
63' 

Adding Eqs. ©-(H), we obtain 
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Thus, the a 7 ln(l/a) term in Eq. (|3|) has the same sign as the a 7 ln 2 (l/a;) one and amounts 
to about 1/3 of the latter. It reduces Az/ by 323 kHz, while the a 7 ln 2 (l/a) term reduces 
Au by 918 kHz. For comparison, we recall that, in the counterpart of Eq. (^) appropriate 
for the muonium HFS, the coefficient of the a 7 ln 2 (l/a) correction reads —8/3 and C = 



281/180 - (8/3) In 2 « -0.2873 Our final prediction for the HFS reads 

Az/ th = 203.39169(41) GHz. 



(9) 
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Here, the uncertainty due to the unknown nonlogarithmic 0(a 7 ) term in Eq. (|3]) is es- 
timated by using the value D = 16.2337T 2 of the analogous coefficient in the case of the 



HFS of muonium [21 ,26] . 



The unknown nonlogarithmic 0(a 7 ) term in Eq. (|3]) receives contributions from three- 
loop QED diagrams with a considerable number of external lines, which are still beyond 
the reach of presently available computational techniques. In this sense, we expect Eq. (|9]) 
to remain the best prediction for the foreseeable future. 

The new theoretical value in Eq. (|9]) exceeds the experimental values in Eqs. (|l|) and 
(@) by approximately 2.6 and 3.5 experimental standard deviations, respectively. Thus, 
the discrepancy between experiment and theory is somewhat reduced by the inclusion of 
the a 7 ln(l/a) term, but is still remains sizeable. 

We may speculate about the magnitude of the coefficient D in Eq. (||). Note that 
two powers of a in the nonlogarithmic 0(a 7 ) term can be of nonrelativistic origin. Each 
of them should be accompanied by the characteristic factor ir, which happens for the 
logarithmic terms. Thus, a plausible estimate of the coefficient D is a few units times ir 2 . 
In order to bring the theoretical estimate into agreement with Eqs. ([I]) and (|2|), we need 
D fh — 1007T 2 and D fa —70ir 2 , respectively. On the other hand, the logarithmic terms of 
the positronium HFS exhibit a structure similar to the muonium case, so that it is not 
unreasonable to expect the nonlogarithmic terms of the positronium and muonium HFS's 
to be of the same magnitude. This would imply a significant contradiction between the 
current experimental measurements and the theoretical prediction. However, although 
this may seem unlikely, one cannot completely exclude the possibility that the residual 
discrepancy will be removed by the inclusion of the nonlogarithmic 0(a 7 ) term. We 
conclude that both the calculation of the coefficient D as well as improved experimental 
measurements are necessary in order to establish or remove the residual discrepancy. Al- 
though there is no conceptual problem on the theoretical side, from the technical point of 
view, such a calculation represents a challenge for QED bound-state perturbation theory. 

Finally we would like to note that the technique developed in Ref. [|13| and here 
can be applied to the analysis of QCD heavy quark-antiquark bound states, where the 
logarithmically enhanced corrections are known to be essential [|T^,^]. 
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